\(\int (d+e x)^2 (a d e+(c d^2+a e^2) x+c d e x^2)^p \, dx\) [2093]

Optimal result

Integrand size = 35, antiderivative size = 95 \[ \int (d+e x)^2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, dx=-\frac {(a e+c d x) (d+e x)^3 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \operatorname {Hypergeometric2F1}\left (1,2 (2+p),4+p,\frac {c d (d+e x)}{c d^2-a e^2}\right )}{\left (c d^2-a e^2\right ) (3+p)} \]

[Out]

-(c*d*x+a*e)*(e*x+d)^3*(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p*hypergeom([1, 4+2*p],[4+p],c*d*(e*x+d)/(-a*e^2+c*d^
2))/(-a*e^2+c*d^2)/(3+p)

Rubi [A] (verified)

Time = 0.05 (sec) , antiderivative size = 124, normalized size of antiderivative = 1.31, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.086, Rules used = {691, 72, 71} \[ \int (d+e x)^2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, dx=\frac {\left (c d^2-a e^2\right )^2 (a e+c d x) \left (\frac {c d (d+e x)}{c d^2-a e^2}\right )^{-p} \left (x \left (a e^2+c d^2\right )+a d e+c d e x^2\right )^p \operatorname {Hypergeometric2F1}\left (-p-2,p+1,p+2,-\frac {e (a e+c d x)}{c d^2-a e^2}\right )}{c^3 d^3 (p+1)} \]

[In]

Int[(d + e*x)^2*(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2)^p,x]

[Out]

((c*d^2 - a*e^2)^2*(a*e + c*d*x)*(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2)^p*Hypergeometric2F1[-2 - p, 1 + p, 2
+ p, -((e*(a*e + c*d*x))/(c*d^2 - a*e^2))])/(c^3*d^3*(1 + p)*((c*d*(d + e*x))/(c*d^2 - a*e^2))^p)

Rule 71

Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Simp[((a + b*x)^(m + 1)/(b*(m + 1)*(b/(b*c
- a*d))^n))*Hypergeometric2F1[-n, m + 1, m + 2, (-d)*((a + b*x)/(b*c - a*d))], x] /; FreeQ[{a, b, c, d, m, n}
, x] && NeQ[b*c - a*d, 0] &&  !IntegerQ[m] &&  !IntegerQ[n] && GtQ[b/(b*c - a*d), 0] && (RationalQ[m] ||  !(Ra
tionalQ[n] && GtQ[-d/(b*c - a*d), 0]))

Rule 72

Int[((a_) + (b_.)*(x_))^(m_)*((c_) + (d_.)*(x_))^(n_), x_Symbol] :> Dist[(c + d*x)^FracPart[n]/((b/(b*c - a*d)
)^IntPart[n]*(b*((c + d*x)/(b*c - a*d)))^FracPart[n]), Int[(a + b*x)^m*Simp[b*(c/(b*c - a*d)) + b*d*(x/(b*c -
a*d)), x]^n, x], x] /; FreeQ[{a, b, c, d, m, n}, x] && NeQ[b*c - a*d, 0] &&  !IntegerQ[m] &&  !IntegerQ[n] &&
(RationalQ[m] ||  !SimplerQ[n + 1, m + 1])

Rule 691

Int[((d_) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Dist[d^m*((a + b*x + c*x^2
)^FracPart[p]/((1 + e*(x/d))^FracPart[p]*(a/d + (c*x)/e)^FracPart[p])), Int[(1 + e*(x/d))^(m + p)*(a/d + (c/e)
*x)^p, x], x] /; FreeQ[{a, b, c, d, e, m}, x] && NeQ[b^2 - 4*a*c, 0] && EqQ[c*d^2 - b*d*e + a*e^2, 0] &&  !Int
egerQ[p] && (IntegerQ[m] || GtQ[d, 0]) &&  !(IGtQ[m, 0] && (IntegerQ[3*p] || IntegerQ[4*p]))

Rubi steps \begin{align*} \text {integral}& = \left (d^2 (a e+c d x)^{-p} \left (1+\frac {e x}{d}\right )^{-p} \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p\right ) \int (a e+c d x)^p \left (1+\frac {e x}{d}\right )^{2+p} \, dx \\ & = \frac {\left (\left (c d-\frac {a e^2}{d}\right )^2 (a e+c d x)^{-p} \left (\frac {c d \left (1+\frac {e x}{d}\right )}{c d-\frac {a e^2}{d}}\right )^{-p} \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p\right ) \int (a e+c d x)^p \left (\frac {c d^2}{c d^2-a e^2}+\frac {c d e x}{c d^2-a e^2}\right )^{2+p} \, dx}{c^2} \\ & = \frac {\left (c d^2-a e^2\right )^2 (a e+c d x) \left (\frac {c d (d+e x)}{c d^2-a e^2}\right )^{-p} \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, _2F_1\left (-2-p,1+p;2+p;-\frac {e (a e+c d x)}{c d^2-a e^2}\right )}{c^3 d^3 (1+p)} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.06 (sec) , antiderivative size = 112, normalized size of antiderivative = 1.18 \[ \int (d+e x)^2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, dx=\frac {\left (c d^2-a e^2\right )^2 (a e+c d x) \left (\frac {c d (d+e x)}{c d^2-a e^2}\right )^{-p} ((a e+c d x) (d+e x))^p \operatorname {Hypergeometric2F1}\left (-2-p,1+p,2+p,\frac {e (a e+c d x)}{-c d^2+a e^2}\right )}{c^3 d^3 (1+p)} \]

[In]

Integrate[(d + e*x)^2*(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2)^p,x]

[Out]

((c*d^2 - a*e^2)^2*(a*e + c*d*x)*((a*e + c*d*x)*(d + e*x))^p*Hypergeometric2F1[-2 - p, 1 + p, 2 + p, (e*(a*e +
c*d*x))/(-(c*d^2) + a*e^2)])/(c^3*d^3*(1 + p)*((c*d*(d + e*x))/(c*d^2 - a*e^2))^p)

Maple [F]

\[\int \left (e x +d \right )^{2} {\left (a d e +\left (e^{2} a +c \,d^{2}\right ) x +c d e \,x^{2}\right )}^{p}d x\]

[In]

int((e*x+d)^2*(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p,x)

[Out]

int((e*x+d)^2*(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p,x)

Fricas [F]

\[ \int (d+e x)^2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, dx=\int { {\left (e x + d\right )}^{2} {\left (c d e x^{2} + a d e + {\left (c d^{2} + a e^{2}\right )} x\right )}^{p} \,d x } \]

[In]

integrate((e*x+d)^2*(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p,x, algorithm="fricas")

[Out]

integral((e^2*x^2 + 2*d*e*x + d^2)*(c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)^p, x)

Sympy [F]

\[ \int (d+e x)^2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, dx=\int \left (\left (d + e x\right ) \left (a e + c d x\right )\right )^{p} \left (d + e x\right )^{2}\, dx \]

[In]

integrate((e*x+d)**2*(a*d*e+(a*e**2+c*d**2)*x+c*d*e*x**2)**p,x)

[Out]

Integral(((d + e*x)*(a*e + c*d*x))**p*(d + e*x)**2, x)

Maxima [F]

\[ \int (d+e x)^2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, dx=\int { {\left (e x + d\right )}^{2} {\left (c d e x^{2} + a d e + {\left (c d^{2} + a e^{2}\right )} x\right )}^{p} \,d x } \]

[In]

integrate((e*x+d)^2*(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p,x, algorithm="maxima")

[Out]

integrate((e*x + d)^2*(c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)^p, x)

Giac [F]

\[ \int (d+e x)^2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, dx=\int { {\left (e x + d\right )}^{2} {\left (c d e x^{2} + a d e + {\left (c d^{2} + a e^{2}\right )} x\right )}^{p} \,d x } \]

[In]

integrate((e*x+d)^2*(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^p,x, algorithm="giac")

[Out]

integrate((e*x + d)^2*(c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)^p, x)

Mupad [F(-1)]

Timed out. \[ \int (d+e x)^2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^p \, dx=\int {\left (d+e\,x\right )}^2\,{\left (c\,d\,e\,x^2+\left (c\,d^2+a\,e^2\right )\,x+a\,d\,e\right )}^p \,d x \]

[In]

int((d + e*x)^2*(x*(a*e^2 + c*d^2) + a*d*e + c*d*e*x^2)^p,x)

[Out]

int((d + e*x)^2*(x*(a*e^2 + c*d^2) + a*d*e + c*d*e*x^2)^p, x)